Europhysics Letters 



PREPRINT 



Extreme events driven glassy behaviour in granular media 

G. D'Anna 1 , P. Mayor 1 , G. Gremaud 1 , A. Barrat 2 and V. Loreto 3 

1 Institut de Physique de la Matiere Complexe, Faculte des Sciences de Base, Ecole 
Polytechnique Federate de Lausanne, CH-1015 Lausanne, Switzerland 

2 Laboratoire de Physique Theorique Unite Mixte de Recherche UMR 8627, Bdtiment 
210, Universite de Paris-Sud, 91405 Orsay Cedex, France 

3 Universita degli Studi di Roma "La Sapienza", Dipartimento di Fisica, P.le A. Moro 5, 
00185 Rome, Italy and INFM, Center for Statistical Mechanics and Complexity, Rome, 
Italy 



Abstract. - Motivated by recent experiments on the approach to jamming of a weakly 
forced granular medium using an immersed torsion oscillator [Nature 413, 407 (2001)], we 
propose a simple model which relates the microscopic dynamics to macroscopic rearrangements 
and accounts for the following experimental facts: (1) the control parameter is the spatial 
amplitude of the perturbation and not its reduced peak acceleration; (2) a Vogel-Fulcher- 
Tammann-like form for the relaxation time. The model draws a parallel between macroscopic 
rearrangements in the system and extreme events whose probability of occurrence (and thus 
the typical relaxation time) is estimated using extreme-value statistics. The range of validity 
of this description in terms of the control parameter is discussed as well as the existence of 
other regimes. 



Introduction. - Granular matter physics |l|,|^] is a very interesting laboratory for address- 
ing open problems of non-equilibrium statistical mechanics, such as the nature of slow glassy 
dynamics Q and jamming 0], the mechanisms of pattern formation, the physics of avalanche 
phenomena. In this context, recent experimental results have shown [^|, the intriguing 
analogy between the way a perturbed granular medium progresses towards complete rest by 
decreasing the amplitude of external taps, and the vitrification of glass-forming materials. 

This slow glass-like granular behavior raises one of the main questions in the framework 
of granular matter, namely the link between the macroscopic response to an external per- 
turbation and the microscopic (or mesoscopic) processes occurring at the scale of the single 
particle. Recently, it has been conjectured |Q that the occurrence probability of rare events 
might explain the origin of the slow dynamic behavior of granular matter. Such rare events 
arise when the largest externally-induced vibration in the granular medium overcomes some 
suitable threshold, triggering in this way a macroscopic "fracture" of the granular solid, that 
is a grain rearrangement. In this paper we further develop this idea and provide a model close 
to the intuitive picture given by experiments and which may be of interest for similar systems 
such as glasses or disordered systems P,p|,pXj[. 
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Experimental results. - We first recall the basic experimental results, and Fig. [l| which is 
in part adapted from previously published work j5| , is used to introduce the method and sum- 
marize the basic data. The dynamic behavior of the perturbed granular medium is deduced 
from noise measurements. The granular noise, denoted |#(/)| 2 , is obtained by detecting the 
irregular motion of a torsion oscillator (Fig. [j] a), deeply immersed into the perturbed granu- 
lar material, and taking the squared amplitude of the Fourier transform of the observed data 
series. The granular medium is perturbed by shaking the container by well isolated taps, or 
by continuous vibrations. An accelcrometer on the container measures the intensity of the 
perturbation, quantified by the reduced peak acceleration T = a s uj^/g, written in terms of 
the amplitude a s and the frequency f s — lu s /2tt of a sinusoidal shaking, and the acceleration 
of gravity g. Below the ideal fluidization limit at = 1, and at low-frequency (below the 
natural frequency of the oscillator), we observe l// 2 noise spectra. 




Fig. 1 - (a): Typical time series 8(t) obtained at T = 0.056 and f s = 113 Hz.. Inset: Sketch 
of the immersed torsion oscillator, (b): Relaxation time t x vs. T for various forcing frequencies 
fs. Dashed line: VFT fit t x = Aexp[B/(T - To) p ] with A = 6.5 x 10~ 4 , B = 0.0072, p = 0.95, 
To = 1.37 x 10 — 4 rri 1 ' 2 . Inset: Diffusion coefficient (arbitrary units) D vs. 1/vT, obtained directly 
from time series such as in (a). 

This low-frequency noise is the structural, or configurational noise. This is easily under- 
stood considering tapping experiments, where after each tap the granular system completely 
stops in a static configuration, which in turn determines an angular position of the immersed 
oscillator. A series of such taps drives the granular medium from one static (or jammed) 
configuration to another, and the corresponding series of static (or " low- frequency" ) angular 
positions 0(t) of the immersed oscillator probes this process, which displays l// 2 noise. Before 
going further we emphasize some experimental aspects: first, when continuous vibrations are 
used, the low- frequency spectrum is the same as the one seen in tapping experiments || . This 
gives us a very rapid and efficient method to measure the noise, and most of the data are 
obtained in this way. Second, our measurements are performed after the sample has been pre- 
pared in a reproducible state. In particular, after having been poured into the container, the 
granular medium is subject to very strong vibrations (large T) whose intensity is progressively 
reduced until reaching the value of T where the measurements are performed. This discards 
the strong compaction effects occurring during the first run on a loose granular system, should 
we have initiated the measurement directly at a low T after pouring. The procedure in fact 
brings the system directly into an almost stationary state which should correspond to the 
reversible curve described in ||. 
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In order to have a deeper insight on the fluctuations of the rotation angles, we have directly 
used the time series of 6{t) (one example is shown in Fig. |l|a) to compute the structure 
factor S(t) — ((6(t + r) — d(t)) 2 ) for different values of T. The obtained linear dependence 
with r allows us to determine a diffusion coefficient D(T), which features an exp(— B'/y/T) 
dependence at low T, as displayed in the inset of Fig. [l]b. Moreover, the l// 2 noise resulting 
from the time series also gives important hints on the nature of the slow glass-like granular 
dynamics: since |#(/)| 2 is proportional to a diffusion coefficient, its inverse at a given frequency 
can be considered to be proportional to an intrinsic configurational relaxation time, r x , i.e., 
l/l^(/)| 2 = C- T x- ln order to obtain the constant of proportionality C, the noise data are 
compared with susceptibility measurements 0, since a peak in the loss factor (i.e. a peak 
in the tangent of the argument of the complex susceptibility) arises when lo p t x = 1, where 
f p = uj p /2tt is the frequency at which the susceptibility is measured. This gives approximately 
C = 27T1500 mrad^s. 

The following important points can be deduced from the experimental measurements. 
First, previous experiments p] provide evidence for a parameter of the form T = b\ZT/ui s = 
by a s /g, with b a constant. This key control parameter determines the "level" of the low- 
frequency noise, as shown in 0]. In other words, whatever amplitude a s and frequency f s of 
the perturbation is, provided T < 1, the noise only depends on T. (In Fig. [j]b we set b = g 1 ! 2 
and the control parameter is T = y^al. Notice that, since we have not changed g, experiments 
only prove that the control parameter is proportional to the square root of the perturbation 
amplitude, i.e., T oc */a7). Another indication for a control parameter like T comes from the 
direct measurement of the diffusion coefficient (inset of Fig. |l|b) which shows an exp(— B' /T) 
form, in agreement with the measurements of t x . The physical signification of this empirical 
control parameter will be discussed in detail below. 

Summarizing, from noise measurements as a function of Y at different shaking frequen- 
cies f s , and having determined the relationship between |#(/)| 2 and t x from susceptibility 
measurements, we obtain the curves shown in Fig. [l]b, namely a plot of the configurational 
relaxation time t x as a function of the control parameter T. Secondly, the data for small T 
in Fig. [j]b have been fitted with the expression t x = Aexp[B/(T — To) p ], which gives in par- 
ticular p — 0.95. Since the observed exponent is close to p — 1, we assume from now on that 
the configuration relaxation time is well described by the standard Vogel-Fulcher-Tammann 
(VFT) expression 

r^-^exp^T-To)]. (1) 

In this VFT form, T is the empirical control parameter playing the role of temperature (which 
does not mean that it is a temperature!), and To is the value of the control parameter where 
the configuration relaxation time scale diverges. (*) 

Model. Before going further, it is important to identify the granular regime we are 
addressing. It is evident from Fig. [l]b that T appears as the relevant control parameter only 
for external shaking with T < 1. In fact, the scaling described by eq.(|lj) ceases to hold for values 
of T which correspond, for each / s , to values of T around unity. This is an important point 
and the model is supposed to be valid for T < 1. In this regime no fluidization occurs and the 
granular medium can be considered in a quasi-solid phase where geometric frustration plays a 



( x ) Actually the data for t x can also be fitted by the simple exponential form t x ~ Acxp[B/T], as well as by 
the form ylexp 



B/ 



f2 



X2 
" L 



, at least for T ^> Tq; we shall come back to this point on the discussion 



of the model. (An experimental estimate of To gives 1.37 X 10~ 4 m 1 / 2 , which corresponds to a displacement 
of about 19 nm. Notice that To is very sensitive to the exponent p, and should be considered only as an order 
of magnitude.) 
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major role. During a periodic vertical shaking of amplitude a s and frequency / s , with T < 1, 
the granular medium moves as a whole, following the movements imposed to the container. 
At the same time the torsion oscillator immersed in the medium is fixed in the laboratory 
reference frame, i.e., it moves with respect to the medium with the same periodicity imposed 
by the shaking. It turns out that it is easier to explain what happens during the shaking 
if one inverts the point of view and considers the granular medium as globally immobile in 
the laboratory reference frame, while the torsion oscillator is subject to a vertical oscillatory 
forcing. Thus, at each period the torsion oscillator penetrates the granular medium for a 
length a s and returns to the original position. In other words, one externally imposes a given 
displacement (or deformation) to the medium, while the stress exerted is not fixed. 

What happens at the microscopic level? During the excursion of length a s , the torsion 
oscillator will try to find its way inside the granular medium by advancing and possibly 
rotating of a certain angle in order to adapt to the reaction of the granular system. We have 
the following picture. For small excursions a s (to be defined below) the medium stays in 
a given static (or jammed) configuration, and it is deformed either elastically or plastically. 
Eventually the medium fractures, and jumps to another jammed configuration. Our immersed 
oscillator detects essentially each fracture-event, and a time-series represents a sequence of 
different jammed configurations. In contrast to an ordinary solid, for which a fracture is a 
unique fatal event, the granular solid is able to reestablish itself in a new jammed configuration 
and to support successive fractures. 

Thus we assume that the penetration/deformation process involves three distinct regimes: 
1) An elastic regime for excursions a s below a limit denoted a°. For a s < a° the system is 
able to absorb the imposed displacement elastically. (The experimental value of a° is material 
dependent and, as discussed below, we identify it with the displacement for which we observe 
the divergence at To, i.e., a° w 10~ 8 m.) In this regime the granular system responds as 
an elastic medium and the torsion oscillator returns to the original angular position at the 
end of the cycle; 2) A plastic regime, in which the imposed displacement is absorbed by the 
granular system by reorganizing the internal stress network, inducing irreversible rotations of 
the angular oscillator position, but without leaving the actual jammed configuration; 3) A 
fracture process in which the system is unable to further absorb the imposed displacement, 
and a macroscopic grain rearrangement (or internal avalanche) is required, resulting in a large 
jump of the torsion angle 0(t). 

We now formalize this three regime process. In the plastic phase the torsion angle under- 
goes sudden (but very small) changes due to the reorganization of the internal stress distribu- 
tion. In a first approximation one can imagine that the variations of the torsion angle could 
be described by a discrete random walk, i.e. one has a gaussian distribution for the torsion 
angles whose variance has to be computed by estimating the number of random walk steps the 
torsion oscillator undergoes for an excursion of length a s . We shall see in the following how 
the hypothesis of a gaussian distribution, partially supported by the results for the diffusion 
coefficient, is not crucial and can be relaxed without changing the main results. 

Since the relaxation law ([l]) implies that below a certain value of a s the relaxation time 
is always infinite, we can deduce the existence of an elastic threshold for a s , i.e. <z°. This 
means that for a s < a® no permanent deformation is produced in the system and the typical 
time for a macroscopic rearrangement is infinite, hence the existence of an elastic phase. For 
a s > a° the reaction of the medium on the torsion oscillator induces irreversible rotations, 
hence the existence of a plastic phase. For these reasons a° can be considered as the length 
of the elementary step of the random walk. 

Now, provided that a s > a°, the number of random walk steps n performed by the torsion 
oscillator for a displacement a s will be proportional to a s — a®. If the elementary angle of 
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rotation at each step is of the order of 9 e i we expect a gaussian distribution for the rotation 
angles with variance a oc 9 e i\fn oc 9 e i \J (a s — a®) /a®. This is a first indication that a parameter 
proportional to ^/a^ oc T (or y/(a s — a°)/ a s K \/T 2 ) cou ld pl a Y the role of a control 

parameter, i.e. a role similar to a temperature in the sense that it determines the variance of 
the stress fluctuations induced indirectly by the advancing torsion oscillator. 

Let us now try to explain the VFT behavior of the relaxation time. This time represents the 
typical time for the system to undergo a macroscopic grain rearrangement. The question can 
be then rephrased as follows: what is the probability that the torsion oscillator, in its angular 
random walk with variance proportional to 6 e i\J (a s — a 2)/°sj w ^ produce a macroscopic 
configurational change such that the system jumps to another static (jammed) configuration? 
Our hypothesis is that the fracture event takes places as an extreme event whose probability 



can be computed in the framework of the extreme order statistics 11 . Fluctuations in the 
torsion angle correspond to a stress redistribution and one expects larger angles to produce 
larger compression of the grain-chains involved. It is natural to expect that there will be 
some threshold value 9 / for the fluctuating torsion angle above which the contact network will 
break and yield to a macroscopic rearrangement leading to a large jump for 9(t). We have 
then to look for the probability that, among the fluctuations of the torsion angle excited at 
each shaking (i.e., at each penetration over a s ), the largest fluctuation will be larger than a 
given threshold 9f. 

Wc define Z n = max(Ai, X2, A„), with Xi — 9ija. We search the probability that 
Pr[Z n > xt] = S n (xf), where Xf = 9f/o~ is the normalized angular threshold. The (cumula- 
tive) probability distribution is Pr[Z n < x] = H n (x), and S n (x) = 1 — H n (x). According to 
standard text-books, for a Normal parent probability density distribution, using the sequences 
given by c„ = V 2 In n — (In Inn + ln47r) / (2-\/2 Inn) and d n = 1 / V 2 In n |ll]], the probability 
distribution H n (c n + d n x) tends, as n increases, to a Gumbel probability distribution |12|j , 
i.e. lim H n (c n + d n x) = exp[— exp(— x)). For large n the probability S n {xf) is given by 

n — >oo 

S n (x f ) = 1 - exp{-exp[-(x/ - c n )/d n }}. 

For (xf — c n )/d n ^ 1, the probability S n {xf) can be approximated by an exponential 
function, i.e., S„(xf) w exp[— (x/ — c n )/d n ]. Using the expressions for c n and d n given 

above, the above condition is verified if n 2 exp n ^ V 47r In n. In our case a oc 

9 e i\/{a s - aP s )/a° s and the approximation is correct if 9 f /9 el > \/2nlnn- %/ " 1 ^ 2 ^ ln ", which 
is a reasonable assumption given the experimental values of the parameters. Using the previous 
approximation the probability for a rare fracture event (corresponding to a configurational 
rearrangement) is 

S n (9 f /a) » n exp \d*L y/^kT, 
V 47r In n 



(2) 



We expect that 9f will depend on the specific tribological properties of the grains, as well as 
on some geometrical properties of the system, such as the grain shape and size distribution. 
Thus, for a given granular system (i.e., fixed n and 9f) the fracture probability is determined 
only by a. The inverse of the probability S n (9f /a) determines the characteristic time for grain 
configuration rearrangements, i.e. the characteristic time of the macroscopic dynamics. One 
has then S n oc t" 1 . Recalling that a is directly related to the empirical control parameter 
T oc y/ch, one finds 



A exp 



BU T 2 - T 



(3) 



which, given the extremely small experimental value of a°, is indistinguishable from (pi) and 
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also very close to an Arrhenius behavior. As already stressed the experimental data do not 
allow to definitely discriminate between the three behaviors. 

It is important to notice that in the present model the gaussian distribution of the torsion 
angle in the plastic regime follows from the random walk analysis. However the precise form 
of the torsion angle distribution is irrelevant as long as it decays faster than any power-law 
(see fl(|). For instance for stretched exponential distributions as exp(— a\6\P) one would 
obtain again a Gumbel distribution for the largest fluctuation where the parameters c n and 
d n would be given, to logarithmic accuracy, by c„ ~ (21nn/a:) and d n = l/c n . This leads 
to the conclusion that, independently of the precise form of the torsion angle distribution (i.e. 
under very mild assumptions on it) one gets the same result for the relaxation time provided 
the variance of the distribution is proportional to y/ (a s — a°)/a°. 

Conclusions. We have proposed a microscopic model which is able to explain some 
experimental results for the relaxation dynamics of a granular medium described The 
crucial hypothesis is that at low T (T < T / where ~ 1 in our experiments) the macroscopic 
dynamics is controlled by extreme events of the stress fluctuations in the system: if the 
variance of these stress fluctuations, driven by the imposed penetration of the oscillator over 
the distance a s , is proportional to yj {a s — a®) /a® (as it turns out if the microscopic torsion 
angle fluctuations are Gaussian), then a macroscopic change of 9(t) can occur only if such 
microscopic fluctuations overcome a threshold, i.e. if an extreme fracture event takes place. 
These simple ingredients immediately lead to an activated-like, VFT-like character of the 
relaxation time. The deviation of the relaxation behaviour from the Arrhenius law is explained 
by the existence of an elastic threshold a®, which is material dependent: for a s < a° the system 
can absorb the perturbation without any rearrangement. 

For r > Tf the situation is quite different: at each cycle most of the induced fluctuations 
in the system are large enough for a macroscopic rearrangement to occur, and the dynamics 
is not controlled anymore by rare events. In this fluid-like phase, structural (configurational) 
changes are not rare events. Recently Philippe and Bideau have identified, using a geometry 
very similar to the one of our experiments, a threshold value of the order of 1.2 above which 
the system is able to reach a stationary state with a density relaxation ruled by a stretched 
exponential law [ Oj . From this point of view we speculate that the threshold value for 
fluidization could also mark the boundary between a glassy region (low T) where the relaxation 
is driven by extreme isolated events and a quasi-liquid region (high T) where the system is able 
to reach a stationary state. It is important to stress that the threshold value for fluidization 
can depend on the geometry of the container. In particular we expect the threshold value to 
increase a lot for highly confined geometries where the effect of the boundaries is strong. This 
could be the case for the Chicago experiments S where a narrow and tall container is used. 
If the effective threshold for this experiment was large (for instance T f ~ 3-^4) this could 
explain why they observed a very slow relaxation even for values of T well above 1. 

The present work represents only a first step in the direction of a better link between the 
microscopic dynamics of granular media and the macroscopic response to an external per- 
turbation. From this point of view, various experiments can be thought of to clarify several 
points: the actual functional form for the relaxation time as well as the dependence of the elas- 
tic threshold a° on the material, the dependence of the control parameter on the acceleration 
of gravity g, the exploration of the fluid-like regime (r > Tf). Finally the experimental setup 
described in this paper could allow the measurement of some thermodynamical properties of 
granular media: in particular the existence of effective temperatures could be investigated 
by a suitable combination of susceptibility and noise measurements, in the spirit of recent 
experiments on laponite [H . Such experiments would open a way towards a comparison with 



G. D'Anna, P. Mayor, G. Gremaud, A. Barrat and V. Loreto: Extreme events driven glassy behaviour in granular m: 



recent theoretical predictions on a thermodynamical approach to granular matter 15, 16| 



* * * 
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